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Abstract 

The stacking order degree of freedom in trilayer graphene plays a critical role in determining 
the existence of an electric field tunable band gap. We present spatially-resolved tunneling spec- 
troscopy measurements of dual gated Bernal (ABA) and rhombohedral (ABC) stacked trilayer 
graphene devices. We demonstrate that while ABA trilayer graphene remains metallic, ABC tri- 
layer graphene exhibits a widely tunable band gap as a function of electric field. However, we find 
that charged impurities in the underlying substrate cause substantial spatial fluctuation of the gap 
size. Our work elucidates the microscopic behavior of trilayer graphene and its consequences for 
macroscopic devices. 
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Graphene has great potential to be used in novel electronics applications [l| due to its 
extraordinarily rich physical properties [2- 5] . Many of these applications require inducing 
a sizeable band gap without sacrificing its high intrinsic carrier mobility While this has 
thus far not been achieved in single layer graphene, multilayer graphene has the possibility 
of an electric field tunable band gap. Trilayer graphene exhibits two natural stacking orders, 
Bernal and rhombohedral. The more commonly found Bernal-stacked trilayer graphene is 
not expected to exhibit a significant field tunable band gap due to its mirror symmetry [l], 



aa 



10 



12 



15[, and experimental evidence thus far has supported this 



16 



. Its zero- 



field low energy band structure behaves roughl 



ike a decoupled stack of Bernal-stacked 



y like a 



12 



15] . However, rhombohedrally- 



bilayer graphene and single layer graphene [l|, |2j 
stacked trilayer graphene behaves differently in an electric field owing to its lack of inversion 
symmetry. In this case, the low energy bands mimic a single sheet of Bernal-stacked bilayer 
raphene, and as such, a large field tunable band gap is expected in ABC trilayer graphene {4, 



12 



15 



20| . Prior work has explored the zero-field band structure of this material via 



electrical transport measurements and found a small band gap opened due to many-body 



effects 



21 



optically 



. The magnitude of the band gap in a tunable electric field has also been studied 



22]. However, direct spectroscopic confirmation of the electronic properties of 



ABC trilayer with independent control of the Fermi energy and electric field has thus far 
been lacking. 

In this letter, we present scanning tunneling microscopy (STM) and scanning tunneling 
spectroscopy (STS) measurements of dual gated exfoliated samples of both ABA- and ABC- 
stacked trilayer graphene on Si/Si02 substrates. Our results provide a direct understanding 
of the behavior of the band gap of ABC trilayer graphene in an electric field, and illuminate 
the effects of local disorder on its electronic properties. 

Trilayer graphene flakes were mechanically exfoliated on a Si02 substrate. Cr/Au elec- 
trodes were written using electron beam lithography. The devices were annealed at 350 °C 
for 2 h in a mixture of argon and hydrogen and then at 300 °C for 1 h in air before being 
transferred to the ultrahigh vacuum low-temperature STM for topographic and spectroscopic 
measurements. Figure 1(a) shows a schematic diagram of the measurement set-up used for 
imaging and spectroscopy of the trilayer graphene flakes. The layer number and stacking 
order of each flake was characterized using confocal Raman spectroscopy mapping prior to 
electrode deposition. Representative Raman signals for both stacking orders are shown in 



Fig. QJc). ABA stacked trilayers show a nearly symmetric 2D peak while ABC stacked tri- 



layers show an asymmetric peak 2lH24j|. Only flakes exhibiting homogenous stacking order 
were used to ensure that local STM measurements were made on a known stacking order. 
The insets to Fig. [He) show spatial maps of the FWHM of the 2D peaks for an ABA- and 
ABC-stacked sample, indicating that the flakes have a uniform stacking order. 

All the STM measurements were performed in ultrahigh vacuum at a temperature of 
4.5 K. dl/dV measurements were acquired by turning off the feedback circuit and adding a 
small (5 mV) a.c. voltage at 563 Hz to the sample voltage. The current was measured by 
lock-in detection. Topographically, we measure a triangular lattice in both types of trilayer 
graphene samples, as illustrated for ABC trilayer graphene in Fig. QJd), which is consistent 
with the expected topography for a multilayer graphene sample. 
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FIG. 1. (color online) Density of states of ABA and ABC trilayer graphene at varying gate voltage, 
(a) Schematic of the measurement setup showing the STM tip and an optical microscope image 
of one of the measured samples, (b) Band structure for ABA and ABC trilayer graphene with no 
electric field and a moderately sized electric field, (c) Raman spectroscopy of both stacking orders 
of trilayer graphene. Left upper inset: Map of the FWHM of the 2D peak for an ABA-stacked 
flake. Right upper inset: Map of the FWHM of the 2D peak for an ABC-stacked flake, (d) STM 
topography image of ABC trilayer graphene showing the triangular lattice. Scale bar is 2 nm. 
Imaging parameters were sample voltage -200 mV and tunneling current 100 pA. 
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We achieve dual control of the Fermi energy and perpendicular electric field by tuning the 
voltage on the silicon back gate and STM tip. The applied electric field is a superposition of 
four fields: the field due to the voltage applied to the back gate, the field due to the voltage 
applied to the tip, the field due to the intrinsic charged impurities in the SiC>2, and the field 
due to the work function difference between the graphene and the tungsten tip. Application 
of an electric field serves to open a band gap for the ABC-stacked trilayer graphene, as 
shown schematically in Fig. []Jb), while the ABA-stacked trilayer remains metallic. Figs. [2] 
(a) and (b) show dl/dV spectroscopy, which is proportional to the local density of states 
(LDOS), as a function of sample voltage and gate voltage for the ABA- and ABC-stacked 
trilayers respectively. Each curve is the average of 100 measurements over a 20 nm by 20 nm 
region. For both stacking orders, features move more positive in sample voltage as the gate 
voltage becomes more negative, due to the shifting of the Fermi energy. However, we also 
apply a larger perpendicular electric field as the gate voltage becomes more negative. For the 
ABA-stacked trilayer graphene in Fig. EJ^a), we see two peaks separated by roughly 30 meV 
(as well as a dip which remains pinned to zero energy). These peaks, marked by black and 
white arrows, move in parallel with changing gate voltage, which suggests they are relatively 



insensitive to the magnitude of the electric field at the densities we probe. We calculate 



the 



band structure of ABA-stacked trilayer graphene using the low energy Hamiltonian j^, 1(3] 
and find relatively good agreement for the evolution of our features with electric field. Our 
calculation shows that while the single-layer like bands move farther from the Fermi energy 



with increasing field, the bilayer-like bands stay rough 



y evenly spaced and exhibit a fairly 



flat band edge which results in our measured peaks 251 ]. 

For the remainder of this study, we focus on the spectroscopic results for the ABC-stacked 
trilayer graphene, where we open a field tunable band gap. For the ABC-stacked trilayer 
graphene in Fig. Efb), we see a dip in the LDOS which grows in width as the gate tunes 
the center of the dip away from the Fermi energy (which occurs around V g = +25 V). We 
attribute this dip to an electric field induced band gap (gap edges marked by black tick 
marks). We also see a peak just above the conduction band edge which does not appear 
on the valence band side (marked with black arrows). This peak is a signature of the van 
Hove singularity at the band edge, and grows in strength as the field becomes larger and 
the band becomes flatter. This electron-hole asymmetry is present in all samples measured, 
but cannot be explained by the low energy Hamiltonian {4] and is not currently understood. 
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FIG. 2. (color online) Density of states of ABA and ABC trilayer graphene at varying gate 
voltage, (a) Experimental dl/dV curves for ABA trilayer graphene taken in regular intervals of 
back gate voltage between +45 V and -45 V. Black arrows indicate the location of the bilayer-like 
conduction band edge. White arrows indicate the location of the bilayer-like valence band edge, 
(b) Experimental dl/dV curves for ABC trilayer graphene taken in regular intervals of back gate 
voltage between +45 V and -45 V. Black arrows indicate the location of the conduction band van 
Hove singularity. Black tick marks represent the band gap edges. All curves are the average of 100 
measurements taken over a 20 nm by 20 nm range. All curves are offset for visual clarity. 



All samples measured showed considerable spectroscopic variation between locations sep- 
arated on the nanometer scale. To counteract this effect, we recorded dl/dV spectroscopy 
curves every two nanometers over a 50 nm x 50 nm region for a variety of voltages on the 
back gate. We then averaged the curves before analyzing the gap size. By averaging the 
dl/dV spectroscopy over even a small region of the sample, we average over spatial varia- 
tions in the band gap due to the varying local electric field strengths resulting from charged 
impurities in the underlying SiC>2 substrate. To understand the magnitude of this effect, we 
recorded the LDOS at a given gate voltage with higher spatial resolution over roughly the 
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same spatial range. We determined the energy of a given feature for each curve; we chose the 
van Hove singularity for the ABC-stacked trilayer graphene since this feature is the sharpest. 
Fig. |3^a) shows the energy of the van Hove singularity for each spectroscopy curve in a 40 
nm x 40 nm region of the sample with the Fermi energy tuned near the charge neutrality 
point (CNP). Similar to the case of single 26M28] and bilayer graphene 29|, |30j, we notice 



spatially coherent puddles of charge on the sample which have no significant correlation with 
topographic features. Fig. [3](b) shows a histogram of the energy of the van Hove singularity, 
from which we see a roughly Gaussian distribution with a FWHM of 10.7 ± 0.8 meV. Sim- 
ilar measurements on different spots on the same sample, and from different samples, give 
comparable charge variations. We do not find significant changes in the charge variation as 
a function of back gate voltage, which is to be expected given that the low energy bands do 
not change in curvature much with changing energy and electric field. By auto-correlating 
the charge puddle map, we are able to estimate a puddle size of 8 ± 0.2 nm for this sample 
spot, as shown in Fig. E|c). 
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FIG. 3. (color online) Charge fluctuation map for ABC trilayer graphene. (a) Experimental dl/dV 
puddle map for ABC trilayer graphene taken atV g = +30 V. The color scale represents the sample 
voltage of the van Hove singularity of the conduction band. The scale bar is 10 nm. (b) Histogram 
of the sample voltage values from the dl/dV puddle map. The red curve is the experimental data, 
and the blue curve is a best fit Gaussian, (c) Auto-correlation of the dl/dV puddle map. The red 
curve is the experimental data, and the blue curve is a best fit exponential. 



At high induced charge density, puddles of increased charge density give larger electric 
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fields and hence larger band gaps. On the other hand, puddles of reduced charge density 
have lower electric fields and hence smaller band gaps. Given a typical puddle size of less 
than 10 nm, we would expect to have more than 25 puddles in our 50 nm x 50 nm scan 
windows. Consequently, averaging all the curves for a map taken with a high induced charge 
density (large electric field) will likely exhibit a gap size corresponding to the average electric 
field of the map, because spectroscopy curves taken in spots with larger than average charge 
densities are compensated for by curves taken in spots with smaller than average charge 
densities. However, this does not hold true for maps taken near the minimum applied 
electric field (near zero charge density). For both positively and negatively charged puddles, 
the local electric field is larger since more charges of either sign induce a bigger field, so the 
gap size for each individual curve can only be larger than the expected gap size given the 
average electric field of the map. Therefore, by averaging over many puddles, at low charge 
density we expect a larger gap size than anticipated. 

To extract band gap sizes, we adopt a fitting procedure in which we fit a V-shaped 
curve separated by a horizontal line to the region directly surrounding the minimum in the 



LDOS [25|. The gap size is equal to width of the horizontal line. The gray circles in Figs. H] 
(a) and (b) represent the gaps extracted from averaging 625 spectroscopy curves taken at 
each gate voltage. To further explore the potential effects of charge puddles on the gap size, 
we also average our data over a smaller spatial area. The exact size of this area does not 



qualitatively change our interpretation of the data [25( ; therefore we present our results for 
a mesh size of 16 nm x 16 nm, which should contain no more than a few puddles. Within 
each of the nine subsections of a 50 nm x 50 nm map, we average all the curves within the 
subsection then fit a gap to that averaged curve. The colored squares in Fig. E(a) show the 
fitted gap size as a function of gate voltage, where each color corresponds to a different 16 
nm x 16 nm area. This data is in good agreement with the averaged data for the entire map, 
with both showing a gap opening on either side of the CNP at around +25 V on the back 
gate. The back gate serves to both tune the Fermi energy and apply an electric field, so the 
observed behavior is to be expected given that when the gate tunes the Fermi energy to the 
CNP, it also tunes the total electric field near its minimum. We measure an electric field 
induced gap as large as 70 meV at negative back gate voltages and as small as 15 meV near 
the CNP. The inset of Fig. Ufa) plots the energy of the conduction and valence band edges 
as a function of back gate voltage. From this, as well as from Fig. [2]^b), we see that the 



majority of the gap opening is due to movement of the conduction band edge. The valence 
band edge appears roughly pinned to zero sample voltage for back gate voltages far away 
from the CNP. Both band edges move roughly in parallel when the back gate is tuned close 
to the CNP, as the gap is small there and the Fermi energy is no longer pinned near one 
band. If we were able to tune our device to higher electron densities, we would expect to 
see the conduction band edge pinned near the Fermi energy. 




Gate Voltage (V) Density (x1 12 cm" 2 ) 

FIG. 4. (color online) Experimental and theoretical band gap size in ABC trilayer graphene as 
a function of back gate voltage and charge density, (a) Experimental gap size for ABC trilayer 
graphene as a function of back gate voltage. Gray circles represent an average of 625 measurements 
over a 50 nm by 50 nm range. Colored squares represent an average of 64 measurements over a 
16 nm by 16 nm range. All data represented in the colored squares is within the 50 nm by 50 
nm range represented by the gray circles. Statistical error bars are smaller than marker sizes. 
Inset: Sample voltage of the valence (red squares) and conduction (blue squares) band edges as a 
function of back gate voltage. Gray circles and colored squares have the same meaning as in the 
main figure, (b) Experimental gap size for ABC trilayer graphene as a function of charge density. 
All color conventions are identical to (a). The black curve is the theoretically expected gap size as 
a function of charge density. 

We note that our conduction band edge moves at a rate of roughly 1 meV in sample 
voltage for every 1 V on the back gate. We use this to relate the gate voltage to induced 
charge density. Fig. IU(b) shows the gap sizes from Fig. IU(a) re-plotted against charge density, 
where zero density corresponds to a gap centered at zero sample voltage. This should account 
for charge puddle fluctuations, and as a result the gap sizes from different spatial locations 
have less variation. The solid black line represents the theoretically expected band gap size 
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as a function of charge density. This was calculated using the low energy Hamiltonian of 
ABC graphene which accounts for trigonal warping terms and interlayer screening (4]. We 
also add a constant charge density offset due to the work function mismatch between the 
tip and the trilayer graphene. For this device, the global CNP was greater than our range 
of back gate voltages (± 45 V), but the local CNP was around V g = +25 V, so there was at 
least a 20 V offset due to the work function mismatch between the tip and graphene 25]. 
We note that this offset cannot be explained by local charge fluctuations alone, given that 
we see a charge variation corresponding to only about 10 V of back gate. When we account 
for this work function mismatch, we find theoretically that we are unable to close the band 
gap for any back gate voltage. Furthermorejthe minimum theoretical gap is larger than the 



expected zero-field gap of roughly 6 meV 



2l| opened due to electron-electron interactions, 



therefore we safely ignore contributions to the band gap due to many-body effects. We 
see very good agreement in both the slope and magnitude of our band gap at high carrier 
density. At low density, we measure slightly larger than predicted gaps. Again, this is likely 
due to the remaining charge variation within our fine mesh. We see that some gaps are 
larger than predicted by as much as about 10 meV, which is to be expected if the data is 
averaged over a few puddles. We note that in cases where a block of data corresponds to 
one large puddle, such as the light blue dots, the experimentally measured gaps are very 
close to the theoretically expected values over all densities probed. 

We have presented local spectroscopic measurements of both the ABA- and ABC- stacked 
trilayer graphene systems. We find that both interlayer screening and relative work function 
difference between the top gate and graphene play an important role in determining the 
magnitude of the band gap in the ABC trilayer. Finally, we see that local charge fluctuations 
strongly modify the size of the band gap, implying the need to control these if a controlled 
global gap is to be realized. 
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SUPPLEMENTARY MATERIAL 



I. CALCULATION OF ABA BAND STRUCTURE 

We model ABA-stacked trilayer graphene as three coupled honeycomb lattices, each with 
two inequivalent lattice sites. The layers have the standard Bernal-stacking scheme as de- 
scribed in Ref. jl|. We account for intra-layer couplings between nearest neighbors, 70, 
coupling between inequivalent lattice sites lying directly above and below each other, 71, 
nearest-layer coupling between the next closest set of lattice sites, 73 and 74, next-nearest- 
layer coupling between equivalent lattice sites, 72 and 75, and the on-site energy difference 
between different layer lattice sites, A jlj]. We also account for the charge density on each 
layer given by n« induced via top (n t ) and back (n b ) gating leading to potential differences 
between the layers. As such, we adopt the low energy Hamiltonian of Ref. [2] 
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(1) 

with A 1)2 (n) = —a\(n 2 + n 3 — n h )\ and A 2i 3(n) = —a\(n 3 — n h )\. f is defined as f(k x , k y ) = 
e ifc x a /v / 3 + 2e" ifcxao/2%/5 cosA; y ao/2, with a = 2.46 A the length of the in-plane lattice vector. 
We take a = e 2 Co/e n with cq = 3.35 A the interlayer distance and k = 2.3 the dielectric 
screening constant corresponding to graphene layers on Si0 2 . We take 70 = 3.12 eV, 71 = 
0.377 eV, 72 = -0.0206 eV, 73 = 0.29 eV, 74 = 0.12 eV, 75 = 0.025 eV, and A = -0.009 
eV jsj. Our sign convention differs slightly from Ref. (2)] for overall consistency. 

We model our silicon back gate as a parallel plate capacitor capable of inducing a density 
= a g(yg — Vo), where a 9 = 7.19 x 10 10 cm -2 V -1 is determined by the gate capacitance 
with 300 nm of Si0 2 , V g is the applied gate voltage and Vo is the shift of the charge 
neutrality point (CNP) due to intrinsic doping. Similarly, we model our tip as a parallel 
plate capacitor with an effective distance of 1 nm. We also add a constant density offset to 
the tip to account for the work function difference between the tip and trilayer graphene. 
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For ABA graphene and a tungsten tip, we estimate this to be equivalent to -10 V on the back 
gate, as this was the approximate difference between the global and local CNP. Combining 
these effects, we estimate our tip induced density as a function of back gate induced density 
as n t = — 0.07n b — a g V s set , where Kffsct is the effective back gate voltage due to the work 
function mismatch. We use this, as well as the self-consistent screening calculations from 
Ref. jsj, to define the layer densities as n 2 = 0.3(nb + n t ) and = 0.58rib + 0.1n t . 

Fig. [STT a) shows a numerical calculation of the density of states as a function of charge 
density, where darker colors represent lower density of states. A linear shifting of the Fermi 
energy caused by the induced charge from the back gate is not included in this calculation. 
The main features to notice are two spikes surrounding the minimum in the density of states. 
These spikes originate from the van Hove singularities of the bilayer-like bands. They are 
separated by about 25 meV, and their relative separation only weakly depends on the charge 
density. We therefore attribute our two peaks in Fig. 2(b) of the main text to the van Hove 
singularities of the bilayer-like bands. The dark region roughly surrounding -25 meV in 
sample voltage has a non-zero density of states, and thus we expect ABA trilayer graphene 
to remain metallic. 




FIG. SI. Density of states simulations for ABA and ABC trilayer graphene. (a) Calculated density 
of states for ABA trilayer graphene. (b) Calculated density of states for ABC trilayer graphene. 
Darker color represents lower density of states. Both calculations include a work function offset 
due to the tungsten STM tip. 
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II. CALCULATION OF ABC BAND STRUCTURE 



Similarly to the case of ABA-stacked trilayer graphene, we model ABC-stacked trilayer 
graphene as three coupled honeycomb lattices, each with two inequivalent lattice sites. The 
layers have the standard rhombohedral-stacking scheme as described in Ref . . We account 
for the couplings described for the ABA-stacked trilayer appropriate for the ABC-stacked 
trilayer and adopt the low energy Hamiltonian of Ref. (4] 
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where all parameters in Eqn. [2] are defined as in the prior section. 

For our sample presented in the main text, we estimate the work function difference to 
be equivalent to -25 V on the back gate. This is an educated guess, since we were unable to 
directly probe the global CNP. However, prior ABC trilayer samples, where both the global 
and local CNPs were accessible, showed similar work function mismatches. Furthermore, 
changing this offset corresponds to a small horizontal shift of the black curve from Fig. 4(b) 
of the main text, and does not qualitatively change our interpretation of the results. Again, 



we use this along with the self-consistent screening calculations from Ref. (4] to define the 
layer densities as n 2 = 0.2(rib + n t ) and n 3 = 0.6rib + 0.2n t . 

Fig. ISlT b) shows the same density of states calculation as in Fig. ISlT a). but for ABC- 
stacked trilayer graphene. Again, the linear shifting of the Fermi energy due to the back 
gate is not included in the calculation. The dark region roughly surrounding meV in 
sample voltage has exactly zero density of states and represents the band gap. We extract 
the theoretical magnitude of the indirect band e; function of charge density from this 

calculation. 
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III. FITTING OF BAND GAPS 



We fit the band gap using a piecewise denned function given by 



A(x — valence) + B for x < valence 
f(x) = <B for valence < x < conduction 

C(x — conduction) + B for x > conduction 



(3) 



where A, B, C, valence, and conduction are fit parameters. The gap size is the energy 
separation between the valence and conduction parameters. We also add the constraint 
that B must be equal to the minimum value of dl/dV (i.e. the gap minimum). Finally, we 
only fit this function over a small range around the gap to avoid influences from the different 
slopes in dl/dV around the van Hove singularity in the conduction band side. The red curve 
in Fig. IS2T a) shows an example of a dl/dV curve at V g = +45 V averaged over a 50 nm by 
50 nm region. The blue curve is the best fit of the dl/dV curve using our gap fit function. 
Fig IS27 b) shows similar results for V g = -45 V. 
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FIG. S2. Example best fit gap determination curves at different back gate voltages, (a) The red 
curve is the average of 625 spectroscopy curves from the same data set as the main text at V g = 
+45 V. The blue curve is the best fit using our gap function and is only plotted over the fitted 
region of data, (b) Same as (a), but for V g = -45 V. 
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IV. VARYING MESH SIZE 



We do not see qualitative changes in the gap fits as a function of mesh size. Smaller 
mesh sizes give the advantage of less charge variation for each set of curves, but come at 
the expense of fewer curves to average. Conversely, larger mesh sizes offer more curves to 
average but come at the expense of averaging over more charge fluctuations. Fig. [S3] shows 
the same analysis as in Fig. 4(b) of the main text, but using a mesh size of 10 nm by 10 
nm. The gray circles and solid black curve are the same is in Fig. 4(b) of the main text. 
The solid blue dots are the fits for the local gap size as a function of charge density. There 
is a larger spread in gap size at a given charge density for the smaller mesh size, but the 
overall magnitude and trend of the data as a function of charge density remains the same 
as for the larger mesh size. Outlier points are the results of bad spectroscopy curves, which 
get averaged away with larger mesh sizes. 
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FIG. S3. Experimental and theoretical band gap size in ABC trilayer graphene as a function of 
charge density. The gray circles and black curve are the same as in Fig. 4(b) of the main text. 
The blue dots represent the experimental best fit gap size for the same data is in Fig. 4 of the 
main text, but with a mesh size of 10 nm by 10 nm. 
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